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Abstract 

The mathematical properties of a nonlinear parabolic equation arising in the modelling 
of non-newtonian flows are investigated. The peculiarity of this equation is that it may 
degenerate into a hyperbolic equation (in fact a linear advection equation). Depending 
on the initial data, at least two situations can be encountered: the equation may have a 
unique solution in a convenient class, or it may have infinitely many solutions. 

1 Introduction 

Modelling the flow of complex fluids is a very intricate problem which is far from being 
solved up to now. Besides studies which aim at improving phenomenological rheological 
models (purely macroscopic constitutive laws), only a few attempts are made to recover the 
rheological behavior of a complex fluid from elementary physical processes arising in its mi- 
crostructure. 

The mesoscopic model which has been proposed by Hebraud and Lequeux in [3] deals with 
simple shear flows of concentrated suspensions. It is obtained by dividing the material in a 
large number of mesoscopic elements ("blocks") with a given shear stress a {a is a real number; 
it is in fact an extra-diagonal term of the stress tensor in convenient coordinates) and by 
considering the evolution of the probability density p(t, a) which represents the distribution 
of stress in the assembly of blocks. Under various assumptions on the evolution of the stresses 
of the blocks which will be described below, the equation for the probability density p(t, a) 
for a block to be under stress a at time t may be written as: 



Xr\[-<7 C ,<t c ]0) . D(p{t)) 



d t p = -b(t) 8 a p + D(p(t)) dl a p - /wu "°' u ^-' p + 5 (a) on (0; T) x R ; (1.1a) 

1 a 

P>0; (1.1b) 

p(0,a) = Po (a) , (1.1c) 
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where for / € L 1 (R), we denote by 

D(f) = ^[ f(a)da. 

J-O J\a\>a c 

In equation ()l.la)) . Xr\[-o- c .<t c ] denotes the characteristic function of the open set R\ [— a c , a c ] 
and 60 the Dirac delta function on R. Each term arising in the above equation (|l.la|) (HL 
equation in short) has a clear physical interpretation. When a block is sheared, the stress of 
this block evolves with a variation rate b(t) = Go7(t) proportional to the shear rate j(t) (Gq 
is an elasticity constant); in this study, the shear rate j(t), and therefore the function b(t), 
are assumed to be in L 1 2 oc (R + ). When the modulus of the stress overcomes a critical value <r c , 
the block becomes unstable and may relax into a state with zero stress after a characteristic 
relaxation time To. This phenomenon induces a rearrangement of the blocks and is modelled 
through the diffusion term D(p{t)) d% a p. The diffusion coefficient D(p(t)) is assumed to be 
proportional to the amount of stress which has to be redistributed by time unit and the 
positive parameter a is supposed to represent the "mechanical fragility" of the material. 

In all that follows, the parameters a, To and a c are positive, and the initial data po in Ijl.lcJI 
is a given probability density; that is 

Po > 0, poGT^R), / Po = l- (1.2) 

Jr 

We will be looking for solutions p = p(t, a) in C^(LlnLl) such that ap belongs to L™(Ll) to 
the nonlinear parabolic partial differential equation . The subscript a refers to integration 
over R with respect to a, whereas the subscript t refers to time integration on [0, T] for any 
T > 0. Note that the average stress in the material is given by 

T (t) = / ap{t,a)da , (1.3) 
Jn 

and therefore the above condition on ap ensures that the average stress is an essentially 
bounded function of time. 

Actually in practice, the shear rate is not uniform in the flow and in order to better describe the 
coupling between the macroscopic flow and the evolution of the microstructure we introduce 
and study in a second paper |Xj a micro-macro model. In this model p is also a function of 
the macroscopic space variables and the average stress defined by 1)1. 3JI is inserted into the 
macroscopic equation governing the velocity of the macroscopic flow. 

In order to lighten the notation and without loss of generality we assume from now on that 
a c = 1 and To = 1. This amounts to changing the time and stress scales. 

The main difficulties one encounters in the mathematical analysis come from the nonlinearity 
in the diffusion term and also and even more from the fact that the parabolic equation may 
degenerate when the viscosity coefficient D{p) vanishes, and this will be shown to may appear 
only when D{po) = 0. This difficulty is illustrated on a simplified example just below and 
also in Section |S] where we discuss the existence of stationary solutions in the case when the 
shear rate b is a constant. 



2 



Let us first of all look at the following simplified model which already includes the difficulties 
we are going to face to in the study of equation Q1.1JI . We consider the equation : 



d t u = D(u(t))dl a u; 
u(0,a) = 5X]-i,i[( fJ ) > 

where X]-l,l[ is the characteristic function of the interval ] — 1, 1[. The initial condition is 
on purpose chosen in such a way that D(u(t = 0)) = 0. The function u = |x]_i lr(o') is a 
stationary solution to this equation and for this solution D(u(t)) is identically zero. But it 
is not the unique solution to (|1.4|) in C 4 °(L^) n L^(L^). It is indeed possible to construct 
a so-called vanishing viscosity solution for which D(u(t)) > for all t > 0, and there are 
actually infinitely many solutions to this equation. (This statement is obtained as a corollary 
of Lemma 14.21 in Section |I] below.) 

As far as equation (jl.lj) is concerned, we show that, in the case when D{po) = and 6 = 0, 
we may have either a unique or infinitely many solutions, depending on the initial data (see 
Proposition 14.11 in Section HJ . 

On the other hand, we are able to prove the following existence and uniqueness result in the 
non-degenerate case when D(p$) > : 

Theorem 1.1 Let the initial data po satisfy the conditions 

Po G L 1 (R) n L°°(1R) , p >0, Po = l and/ |cr|p <+°o, (1.5) 

Jn Jn 

and assume that 

D( Po ) > . 

Then, for every T > 0, there exists a unique solution p to the system 1)1.1(1 in Lf°(L^ nl^)n 
L 2 t (Hl). Moreover, peLf^D C t °(L^ nL 2 a ), f R p(t, a) da = 1 for all t > 0, D(p) G C t ° and for 
every T > there exists a positive constant v(T) such that 

min D(p(t)) > u(T) . 

0<t<T 

Besides op G Lf(L l a ) so that the average stress r(t) is well-defined by hi. 3(1 in Lf 3 . 

The first step toward the existence proof of solutions to (|l.lj) will consist in the study of 
so-called vanishing viscosity approximations, which are the unique solutions to the family of 
equations 

d tPe = - b(t) d a p £ + (D(p £ (t)) + e) d 2 aa Pe - XR\[-i,i] Pe + D{Pe{t)) 5 (a) ; (1.6a) 
Pe > ; (1.6b) 
p e (0,.) = Po (1.6c) 

(recall that we have rescaled the time and stress units to get To = 1 and a c = 1). Sectional 
below is devoted to the proof of the following 
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Proposition 1.1 (Existence and uniqueness of vanishing viscosity approximations) 

Let T > be given. We assume that the initial data satisfies the same conditions ()1.5j) as in 
the statement of the theorem. Then, for every T > and < e < 1, there exists a unique 
solution Pe to JISI) inLf{Llr\Ll)f\L 2 t (Hl). Moreover, p £ G L^nC^L* nL l ) , D (Pe) G C°, 

p £ = 1 , (1.7) 



< p e < [|po||igp + ^ 



(1.8) 



and /or even/ T > ; i/iere exist positive constants C\{T,po), C2(T,po) and C$(T,po) which 
are independent of e such that 



o<t<T . 
sup 

0<t<T 

and 

r-T 



sup / \a\p £ < Ci(T,p ) , (1-9) 
<t<TJR 

sup / p 2 < C 2 (T, P0 ) , (1-10) 

<t<T JR 

/' (e + D( Pe )) f \d aPe \ 2 < C 3 (T,p ) . (1.11) 



Theorem (jl.ljl is then proved in Section[3]whilc the degenerate case is investigated in Sectional 
Lastly, the description of stationary solutions in the constant shear rate case is carried out in 
Section \5\ 



2 The vanishing viscosity approximation 

This section is devoted to the proof of Proposition 11.11 



We begin with the following : 

Lemma 2.1 (Uniqueness) Let po satisfy (jl.2ft . Then for every T > and < e, there 
exists at most one solution p £ to (|1.6|) in Lj°(Lj. n L 2 ) n L 2 {Hl). Moreover, p £ G C t °(L 2 ) 
(thus, the initial condition makes sense) and 

[ Pe = l, (2.1) 
JR 

for almost every t in [0, T]. 



Proof of Lemma I2.lt We begin with proving that every solution to (|1.6|) in Lf{L\ n L 2 ) n 
L 2 (H^) satisfies (|2.1|l . We fix R > 1 and we consider a cut-off C 2 function 0^ = 4>r{°~) with 
compact support which is equal to 1 when < |o~| < R and to when \a\ > 2R and such 
that 

\<P'r\ < § , (2-2) 
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where here and below C denotes a positive constant that is independent of R. Notice that 4>' 
is equal to on ] — oo, — 2R], on [-R, R] and on [2R, +00 [. 

Now, we multiply IJl.fiaj) by (pR and integrate over [0, t] X R to obtain 

[ p £ {t)4> R - [ Po4>R = - [ b(s) [ daPe(s) 4>R ~ [ (D(p £ (s)) + e) [ d a p e {s) (j>' R 

Jr Jr Jo Jr Jo Jr 



t r 1 ft 



'0 J\<t\>1 a Jo 

We bound from above the terms on the right-hand side as follows. First, we have 



JR 



b( S ) / d aPe {s)<t>R < / \b( S )\ / Pe{s)\4>' l 



J R\ 

'0 JR 

< ~ f \Ks)\ [ p e (s) < 2 

n Jo JR<\a\<2R n 



e\\Ll \WrWl* 



thanks to l(2~2")) and using that p £ G Lf{L\) and be L\. Next, 

I (D(p £ ) + e)\ f d a p £ (j)'\ < (e + a||p e ||x?°(xi)) f PaP, 
Jo Jr Jo 

CVt _C_ 
thanks again to (|2.2jl . Cauchy-Schwarz' inequality and since d a p £ is in L 2 a . Finally, 

0<- f D(p £ ) - [ [ p £ <\>R =11 Pe(l-<I>R) 
a Jo Jo J\ff\>l Jo J\a\>l 

< T / Pe , 



JO J\a\>R 

and the right-hand side goes to as R goes to infinity since p £ is in L^{Lyj. All this together 
yields 

/ Pe(t) = lim / p £ (t)cj)R = lim / p 4>R= / Po = 1 , 
Jr R-++°°Jr R-^+°°Jr Jr 

for almost every t in [0, T]. In particular, this implies that D(p £ ) < a . 

Let us now argue by contradiction by assuming that there exist two solutions p\ and P2 to 
(jl.6a|) corresponding to the same initial data P q. By subtracting the equations satisfied by p\ 
and P2 respectively, we obtain 

8 t q = - b(t) 8 a q + D(q) % ffPl + (D(p 2 ) + e) d 2 aa q - X R\[-l,l] Q + ^ 5q((t) ; (2.3) 
9(0,(7) =0, 

where q = pi — P2- We multiply (J2.3D by q and integrate over R with respect to a to obtain, 
after integrations by parts, 

\t+ I <? + ( D &)+ £ ) I l^l 2 + / I 2 = —Q(t,0)-D(q) [ d aPl d a q. (2.4) 
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We first remark that since pi = / P2 = 1 thanks to ()2.1|) . we get 

Jr Jr 

\ d (q)\ = a \ <?| < »V2\\q\\ L 2 , 

V| CT |<1 

with the help of Cauchy-Schwarz' inequality. Next, using the Sobolev embedding of ff 1 (]R) 
into we bound from above the terms on the right-hand side in the following way : 

\£^-q{t,0)-D(q) [ d aPl d a q\ 
« Jr 

< V2\\q\\ L 2 \\q\\L°° + V2a\\q\\ L 2 / \d a pid a q\ 

Jr 

< 11^11^2 ( ||g|||2 + \\daq\\ 2 L 2)2 + V2 a \\q\\ L 2 ||<9 CT pi|| L 2 \\d a q\\ L 2 

^ 1 II ll2 , a2 |in2 ii a m2 i e n ||2 : ii a n2 



Therefore, comparing with Q2.4I) we deduce 

- 2jt U\ L i < (- + T R^II^ + 2)Nb 



Finally, by applying the Gronwall lemma, we prove that \\q\\ 2 L 2 < 0, thus q = 0. The umque- 
ness of the solution follows. 



Remark 2.1 T/ie same proof shows that if there exists a solution to (jl.lj) in L^°(L^. n L 2 ) Pi 
L 2 (H^) such that info<t<r D(p(t)) > 0, £aen it is unique in this space. 

We now turn to the existence part in the statement of Proposition 11.11 From now on we fix 
a positive constant e < 1. The proof of Proposition 11.11 will be carried out by the Schauder 
fixed point theorem. For given positive constants M(> e) and R, we introduce T> Ej m and Yr 
two closed convex subsets of respectively L 2 and L 2 a as follows: 

V £yM = {a 6 L\\ e <a< M} 

Yr = {p e Ll„, p > , sup I \a\p < R). 

o<t<T Jr 

To simplify notation we denote 

^(x) = ^exp(-^) if > ; 

We first prove the following 

Proposition 2.1 Let T > and Ze£ po G L 2 (R) sitc/i iaa£ po > 0. Then, for every a in T) e .M 
and q in Yr, there exists a unique solution p in L^(L 2 ) n L 2 (H^.) to 

d t p(t,a) = - b{t)d a p{t,a) + a{t)dl a p{t,a) - XB\[-x,x]W)p(f,a) + ^ S (a) ; (2.5a) 
p(0,a) =p (a) . " (2.5b) 
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Moreover, p G C 4 °(L 2 ), p is non-negative and 

P- < P < P+ , (2.6) 

with 

r+oo 

P-{t,a) = e~ t \ po(o-')<p^-(a-o-' - X (t))da' (2.7) 



oo 



and 



/+oo 
fo(cr') <p rr7r-(? - a' - x(t)) da' + 
-oo \J 2 Io a 

i r* 



+- [ D(q(s)) if ,—(a - x(t) + x(s)) ds , (2.i 



a Jo 

where x(t) = fQ°(s)ds. In addition, 

i. If po G L°°(R), then p is in Lf^ and 

o< P < UpoIIl- + 4^t= • ( 2 - 9 ) 

V 71 " v e 

M. // J R |<t| po < +°° f thus po G L X (R) i/ien |<r| p G L^°(L^). More precisely, we have 

r r 2 r 

sup / |cr|p < / kbo + ^11^1^2(0.^ ||p || L i + —T 3/2 ||6|| L 2 (0iT) 
o<t<TJn Jn 6 

+ ^ (MT) l/ 2 | bo|Ul+ i^ r 3/ 2 . (2 . 10) 

V 71 " 3 Vtt 

Moreover p G C^(L\) and D{p) G C t °. 

Proof of Proposition I2.lt Let us first observe that for every q in Yr, -D(g) G since 

< £>(g(t)) < a / |<r| q < a R , (2.11) 
J\a\>l 

for almost every t in [0,T]. Therefore the source term D(q(t)) 5q(o~) in (|2,5a[) is in Lf{H^ 1 ) 
and the existence and the uniqueness of a solution p G C^{I^) n Lf(H^) to the system (|2.5jl 
is well-known (see for example [2])- In particular, the initial condition makes sense. Owing 
to the fact that the source term is non -negative, the proof that p > is also standard (see 
again j5j). 

We now check the pointwise inequality (|2.6j) . 

This is ensured by the maximum principle with observing that p_ and p+ given respectively 
by (|2.7|) and (|2.8|) are the unique solutions to the systems 

d t p- = -bd a p- + adleP- - p_ ; , . 

p_(0,<7) = p (a), 



and 



d t p+ = -h d a p + + a dl a p + + 5 (cr) ; 



P+(0,cr) = po(cr) 



(2.13) 
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respectively. We now turn to the proof of statement i. and assume that po belongs to L°° (R) . 
Then, using the two facts that for every v > 0, J R tp v = 1 and tp v < ()2.9j) is easily 

deduced from p < p+ with the help of (|2,11|) and since a > e. 

Suppose now that J" R \o\ po < +00. This together with the assumption po G L 2 (R), guarantees 
that po G L X (R) (see also below). Using (|2.6j) again, we now have 



<r\P < / W\P+ 
R JR 



+ 



-/ D(q(s))([ \a\(p / -—(a- X (t)+x(s))da)ds 
a Jo Wr \J 2 L a > 



PoW) \a + a' + x(t)\ /TTt-(o-) dado 7 
'r^r v 2 Jo a 

+~ jT* £>(g(a)) ( jf |a + ( X (t) - x(s))\ V f^(?) d °) ds 



1 



< / \<T\po(a)da+\x(t)\\\po\\ L i + - \ X (t)-x(s)\D(q(s))ds 
Jn " Jo 



+ 4=( fa) 1 ' 2 \\p4 L i + * f D(q(s))( f af* ds , (2.14) 



V™ Jo Jo 

since J R \a\ ^Pu( cr ) da = (2/tt) 1 ^ 2 v and J R y? w = 1. With the help of (|2.11|) and observing that 
\x(t) - x(s)\ < y/t-s ||6|| L 2 (0jT) , we then deduce (f2~TU|) . 

We now use this bound to check that p G C°(Lj.) and D(p) G C 4 °. Indeed, for any t, any 
sequence t n in [0,T] which converges to t and A > 1, we have 

P(tn)-P(t)\ = [ \p(t n )-p(t)\+ [ \p(t n )-p(t)\ 
R ' J\cr\<A J\v\>A 



< V2A 

< V2A 



(jjp(t n )-p(t)\ 2 ) +J j R W\(\p(tn)\ + \p(t)\) 

[ \p(t n ) ~ p(t)\ 2 Y /2 + 4 ^p / \a\\p(t)\. (2.15) 
Jn J A o<t<TJK 



For any fixed A the first term in the right-hand side goes to as n goes to infinity since 
p G Cj^L 2 ) and then the second term is arbitrarily small as A goes to infinity. The same 
argument yields the continuity of D(p(t)) with respect to t. 



The following proposition aims at checking the required assumptions to apply the Schauder 
fixed point theorem. 

Proposition 2.2 Let Tf > be given. We assume that 

Po G L*(R) n L°°(R) , p >0, Po = l and/ |cr|p <+°o. (2.16) 

Jn Jn 
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Let < e < 1, R = l + J* R jcr| p and M = 1 + 2 a. We define 



9 

T c = — 
25 



2 VI + 2a 



- 2 



(2.17) 



Then, for every T < min (-^; T c ), the function T : (a;q) i— ► (D(p) +s;p), with p being the 
solution to the system (|2.5|) . maps £> £ ,m x Yr into iise//. Moreover T is continuous and 
T(T> £ ^m x Yr) is relatively compact in L 2 (0, T) x -L 2 0-. 

Proof of Proposition 12.21 

Step 1: T is well-defined. 

According to Proposition 12.11 p is in Cf(L^) and D(p) S Cj . We now prove that with our 
choice for M ( which ensures that e+D(po) < 1+a < M), D{p)+e S P e K . For this, we again 
use the inequality p < p + , the definition (j2.8|) of p+, the rough estimate Jj CT | >1 < Jr V 9 ^ = 1 
and Q2.11JI to obtain 

sup D(p(t)) < sup D(p + {t)) < a + aRT < 2 a , 

0<t<T 0<t<T 

for T < It only remains now to check that sup / \o~\p < R. We thus go back to (f2~T0l) 

0<t<TjR 

and observe that this condition holds provided 

T < max{t > 0; ||6|| L 2 (0iT/) Vt (1 + — t) + -^=- + ^= < 1} . 

Since we already have demanded that t < T < a sufficient condition is then 

^[- 3 W\^,T, ) + ^-\ <1, 

which reduces to T < T c with T c given by (j2,17|) . 



Our next step will consist in establishing a priori bounds on p in Lf°(L 2 ) n L\(H\). 



Step 2: A priori bounds. 

If we multiply equation (|2,5a|) by p and integrate by parts over 1R with respect to a we easily 
obtain 

^/p 2 + a/W<^>p( t ,0). 

Since from the Sobolev embedding of ff 1 (Il) into L°°(Il) and the bound (|2,11|) on D(q) we 
get 

|^p(i,0)| < R\\p\\ L ~ 
a 

< R(M\% + \\dM\%) k 
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we may write 

We recall that a > e and we apply the Gronwall lemma to obtain 

sup \\p\\l 2 < e eT (||p |li 2 + ■ (2.19) 

0<i<T 



We now return to (|2.18j) and integrate it over [0; T] to obtain 

e Rj»||? a < Upollia (1 + eTe £T ) + — (1 + eTe eT ). (2.20) 



Step 3: T/ie function T is continuous. 

We consider a sequence (a n ; g n ) in T> £i m x Yr such that a n converges to a strongly in L\ and 
q n converges to q strongly in L^ a , and we denote T(a n ;q n ) = (D(p n ) + e;p n ). We have to 
prove that p n converges strongly to p in L^ a and D{p n ) converges to D(p) strongly in L|, 
with (D{p) +e;p) =T(a;q). 

In virtue of (|2.19j) and (|2.20() . the sequence p n is bounded in Lf(Ll) n L^(H^). Then, <9 CT p n 
is bounded in Lj°(ff~ 1 ) and d% a p n is bounded in Ll{H~ l ). Since a n d^ a p n is bounded in 
Ll{H~ v ), b € £| and D(q n )5 is bounded in Ll{H~ l ), d t p n is bounded in L^H" 1 ). This 
together with the fact that p n is bounded in Lf(H^) implies that, up to a subsequence, p n 
converges strongly towards p in L^(L^ OC(T ) (the convergence being weak in L^(H^.)) thanks 
to a well-known compactness result [I]. In particular, p ra converges to p almost everywhere. 
Thus p > and by the Fatou's lemma, < R almost everywhere on [0;T]. Hence p 

belongs to Yr. We are going to show that the convergence is actually strong in L^ a . 

In virtue of (|2.9|) in Proposition 12. H we dispose of a uniform a priori bound on p n in 
(hence also on p). For the strong convergence in a we then argue as follows. For any fixed 
positive real number K, we have 

\Pn~P\ 2 < / / \Pn~P\ 2 + / \Pn~P\ 2 

R JO J\a\<K JO J\cr\>K 



f T f ,2 /,, , , , \ 2 RT 

< / / \Pn-p\ + {\\Pn\\L™ + \\P\\L~ 

JO J\a\<K 



K 



owing to the fact that p n and p belong to a bounded subset of Yr n L^° a . We then conclude 
by letting n next K go to infinity. 

We now prove that D{p n ) converges to D{p) strongly in L\. We shall actually prove that 
D{p n ) converges to D(p) strongly in L\ and then use the fact that D(p n ) is bounded in Lf°, 
in virtue of (|2,11|) and because p n lies in Yr. Let us fix K > 1. Then, we have 




D( P )\ = / | / (Pn-P)\ 
Jo J\a\>l 

< f f \Pn~p\ + ^ / / kl (\Pn\ + \p\) 
Jo Jl<\a\<K A Jo J\a\>K 

< / / K - P\ + —r^ , (2.21) 

Jo JK\cr\<K A 
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because p and p n belong to Yr. Since p n converges to p strongly in L\(L\ QCo ), we conclude 
that D(p n ) converges to D(p) in L\ by letting n next K go to infinity in (|2.21j) . 

In order to pass to the limit in the equation satisfied by p n (thereby proving that (D(p) + 
e;p) = T(a; q)), we now observe that the strong convergence of q n to q in a , together with 
the argument in (|2,21|) above shows that D(q n ) converges to D(q) strongly in L\. It is then 
easily proved that p is a weak solution to (|1.6a|) and since p is in Lt(H^) it is the unique 
solution to (|2.5|) corresponding to a and q. In particular, the whole sequence p n converges 
and not only a subsequence. 

Step 4: T(T> £ x Yr) is relatively compact. 

Let (D(p n ) + e ;p n ) = T(a n ; , q n ) be a sequence in T(V £j m x Yr). We have to prove that we 
may extract a subsequence which converges strongly in Li x h\ a . Exactly as for the proof of 
the continuity, the a priori estimates 1)2.19(1 and (|2.2Uj) ensure that the sequence p n is bounded 
in Lf 3 (L^) n Ll(H^). Since \<r\p n is bounded L^{L l a ), we can mimic the argument in Step 3 
above to deduce that up to a subsequence the sequence p n converges to some p in Yr strongly 
in L\ a and that D(p n ) converges to D(p) strongly in . 

We are now in position to conclude the proof of Proposition ll.il 

Let Tf > and < e < 1 being given. We are going to prove the existence of a unique 
solution on [0; Tf]. 

Being given an initial data po which satisfies (|1.5fl , existence of a solution p £ is ensured from 
Proposition 12 . 21 bv applying the Schauder fixed point theorem on "short" time interval [0; T\] 
with T\ = min(-^-, T c ) and where R± = 1+/ R W\pa ■ This solution is uniquely defined in virtue 
of Lemma 12.11 and we know from (|2.1|) that J^p £ (Ti) = 1. Moreover from Proposition 12.11 
Pe(Ti) £ and by construction j* R \a\p £ (Ti) < R\. Therefore p £ {T\) satisfies the same 
conditions (|2.1H|) as p$. Then, repeating the same argument we may build a solution to (jl.fij) 
with initial data p £ (Ti) on [Ti; T2] with T2 = min (-^-, T c ) , where R2 = Ri + 1 = J R \o~\ po + 2. 
Thanks to the uniqueness result (Lemma l2.1|) . if we now glue this solution to p £ at t = T\ we 
obtain the unique solution to (|1.6|) on [0;T\ + T2]. It is now clearly seen that for any integer 

n > 1 we may build a solution to Ifl.fijl on [0; Y^i<k<n Tk] with T& = min [{k+ J R \o\ Po) 1 ; T^j . 
Since J2i<k<n^k obviously goes to +00 together with n, existence (and uniqueness) of the 
solution p £ to (|1.6|) is obtained on every time interval. 

For the proof of (|1.8|) we argue as for the proof of (|2.9|) in Proposition 12.11 Defining pf 
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as in 1)2.8)1 with a replaced by D(p e ) + e and D(q) by D(p e ) we obtain 

< p £ <pf 



Then 



< \\po\\l°° + 

< \\po\\l°° + 

< [|Po||i°° + 

< IIpoI|l°° + 



i r* D(p £ ( s )) 



: ds 



"V^o 2 Je + f t ,D(p e ) 



1 



a y7T 
1 



e + / D( Pe ) - 



a y 7T y 



D(p £ 



7T 



R 



llPelUgP / 

Jn 



Pe 



from which (jl.lUj) follows gathering together 1)2.1)1 and (|1.8[) and, with the notation of the 
proposition, 

C 2 (T,p ) = \\po\\l°° + ^ X 1 



The proof of (jl.9j) follows the same lines as the proof of 1)2.14)1 . Indeed, we again use the 
pointwise inequality p E < pf and replace D(q) by D(p £ ) (< a) and a by D(p £ ) + e (< a + 1) 
in 1)2 .14|) and use (|2.16|) to deduce 

sup / \a\ Pe < [ H P0 + V^(^^ + ||6|| i2(0iT) ) + |T^(l + ^^), (2.22) 
o<t<TjR Jn v v 71 " ' J v V 71 " y 

whence p. 9)1 with Ci(T,po) being the quantity in the right-hand side of (|2.22j) . 

In order to prove 1)1.11)1 . we multiply 1)1. 6a)) by p e , and we integrate by parts over R with 
respect to a to obtain 

~ / p\ + (D(p £ ) + e) [ \d a p £ ? + / p\ = ^p £ (t,0) . (2.23) 

2"tjR JR 7|er[>l a 

We use the L°° bound ()1.8j) to bound the right-hand side and we integrate ()2,23)) with respect 
to t over [0; T] to deduce (jl.ll)) with 

c 3 (r J po) = [|po||icc(i + T)+ ^r 3 / 2 , 

2 v vr 

using that ||j>o|||a < IM|r°° JrPQ- 
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3 The non degenerate case: D(po) > 



The main result of this section corresponds to the statement of Theorem 11,11 and fully de- 
scribes the issue of existence and uniqueness of solutions to the HL equation (|l.lj) in the 
non-degenerate case. It is summarized in the following : 

Proposition 3.1 Let po satisfy (jl.5|) . We assume that D(po) > 0. Then, the HL equation 
(jl.lj) has a unique solution p in C t ° \L%) D L\ (H^) and p is the limit (in L t 2 loc (L^)nC t ° loc (L^)J 
of (p e ) when e goes to where p £ is the vanishing viscosity solution whose existence and 
uniqueness is ensured by Proposition Moreover, p E Lf^ Pi Cf(L^), ap E Lf{L\) and 
J^p = 1. Furthermore, D(p) E C t ° and for every T > there exists a positive constant v{T) 
such that 

min D(p(t)) > v{T) . (3.1) 

0<t<T 



We begin with proving the following : 

Lemma 3.1 We assume that po satisfies 1)1.5)) . Then, if D(po) > 0, D(p e )(t) > for every 
t E [0, T], with p e being the unique solution to ()1.6)) provided by Provosition \l.l\ and, actually, 
for every T > there exists a positive constant v{T) such that 



mm T D(p e (t))>v(T), (3.2) 



for every < e < 1. 



Remark 3.1 iVote i/iai f/izs bound from below is independent of e, but it comes out from the 
proof that it depends on po and on the shear b. 

Proof of Lemma 13. H The proof relies on the bound from below in ()2.6j) that we integrate 
over \a\ > 1 to obtain 

D{p £ (t)) >a I p->ae- 1 f p (a') ( f if , (a-a' - x (t)) da) da' . (3.3) 

J\a\>l JR Wh>1 ^I (D(Pe)+e) J 

Let us define K x = [-1 - x(t), 1 - The function a i-> ip y 2 jt (D(p )+e) ( a ~ a ' ~ 

is a Gaussian probability density with mean a' + and squared width 2 JjJ (D(p £ ) + e). 
Therefore, for every a' E II \ -fT x , we have 



V / o rW v. M ~ °' ~ *(*)) da - \ 
H>i V 2 / (- D ^)+ £ ) 2 



which implies 



(El > ^e" T / Po = %~ T I po 

2 iR\Jf, 2 J| CT+X (t)|>l 

In the zero shear case (6 = 0, thus X = 0) the proof is over and 



min D(p(t)) > -e~ T D(p ) 
o<t<T yFy " ~ 2 y ' 
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In the general case, a strictly positive bound from below is available as long as the support 
of po is not contained in K x . We thus define 

t* = inf \t > 0; / po = j . (3.4) 

1 J\<r+X(t)\>l } 

Then < t* (t* possibly even infinite), the support of po is contained in [— 1 — 1 — x(£*)[, 
and for every T < y, (j3.2|) holds for some positive constant v\(T) defined by 



vi( T ) = T; e T min / po- (3.5) 



0<t<T 



k+x(*)l>i 



It is worth emphasizing that this quantity is independent of e. If t* = +oo, the proof is over 
and v(T) = i>i(T) fits. Let us now examine the case when t* < +oo and T > y. 

We go back to (|3.3j) . take t in [y;T] and denote x = j^{D{p £ ) + e) for shortness. Then 

•i-x(f) 
'-i-x(t') I 



Po(*0 ( / ¥y^( CT - CT ' - *(*)) da) da' 
-i-x(t*) v Jkl>i y 



ae / Po(c') I / ~ — f= — f= do" I dt/ 



-i-x(t*) 



|cr|>i ^\f^\fx 



-7=e / po(ff') / o 7= rf(7 + / o /- dtJ dcr ' 

V 71 " J-1- X (t*) \J-oo 2 V X Jl+a'+x{t) 2 V X 

=e" T / p (a') / e"' (it + / e~* dt da' 



i-x(f) 



~ W-l-v(t*) / \ h~X^)+x(t) I 2+ X (t»)- X (t) 

^ / \ ^2** V 2 ** n(**/2) 

> -^=e~ T min | / e~* 2 di + [ e~ t2 dt\ , (3.6) 

v/7T f/2<t<T \ I 2-x(t*)+x(t) / 2+ x (t*)-x(t) I 

\ V 2 **"i(**/2) V at * "i(**/ 2 ) / 

since Jj^ffjPo = 1 and x > ^ 12 D(p e ) > t*ux(t* /2)/2 thanks to (j33|) . The proof of 
Lemma 13. II then follows by defining 

r/(T) = min( i / 1 (r);z/ 2 (r)) , 

with v\(T) given by (|3.5() and v% (T) being the positive quantity in the right-hand side of 
(|3~6l) . that is 



v 2 {T) = ^=e' 1 min I / e~* 2 di + / e"* 2 dt 

/7T t* /2<t<T \ I 2- x (t*)+x(t) / 2+ x (t*)-xW 

' ^2 i* i/ 1 (i*/2) i/ 1 (i*/2) 

Proof of Proposition 13. H 
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We first go back to the proof of the bound on d a p £ and more precisely we look at Q2.23JI 

and observe that in virtue of (|3.1j) 



v{T) / / \d aPe \ 2 <C 3 (T,p ) . (3.7) 



\"t l'e\ 2 

R 

Let now e n denote any sequence in [0, 1] which goes to as n goes to infinity. To shorten the 
notation we denote hy p n instead oip Sn the corresponding sequence of solutions to (|1.6|) . With 
the above bound l|3.7j) on p n and Ijl.lOj) . we know that p n is bounded in L 2 (H^) independently 
of n. Moreover thanks to (|2.1|l and (jl.8[) p n is bounded in L^°(L^ n and we also dispose 
of a uniform bound on J R \cr\p n in virtue of ()1.9j) . Therefore arguing exactly as in the proof of 
Proposition 12.21 (Step 4) where we have proved that the mapping T is relatively compact in 
L 2 x L 2 a we show that p n converges to some p strongly in L 2 a and D(p n ) converges to D(p) 
in L\. Then p is a solution to the initial problem (|L"Tj) in L 2 t {R\) n Lf{L\ ni™), J R p = 1 
and J* R \cr\p < +oo. Moreover, 

inf D(p(t)) > v(T) . 

0<t<T 

This non-degeneracy condition on the viscosity coefficient ensures that there is at most one 
solution to in L 2 (H^) n Lf°(L 2 ) (this follows by an obvious adaptation of the proof of 
Lemma 12. II to this case). Therefore the limiting function p is uniquely defined and does not 
depend on the sequence s n . Moreover the whole sequence p n converges to this unique limit 
and not only a subsequence. 

As a conclusion of this subsection let us make the following comment which is a byproduct 
of Proposition 13.11 Let p be a solution to (|l.lj) in C t °(Lj. n L%), then as soon as D(p(t)) is 
positive for some time t it remains so afterwards since the solution can be continued in a 
unique way beginning from time t. 

4 The degenerate case : D(po) = 

Throughout this section we assume that D(po) = and therefore the support of po is included 
in [—1; +1]. Assume that we dispose of a solution to (|l.lj) in C^(L l a n L%). We may define 
U £ R+ U {+00} by 



max 



{*>0;^ D(p) = o} . (4.1) 



According to the comment at the end of the previous section for every t > i*, D(p(t)) > 
while D(p(t)) = for all t in [0;t*]. On [0;t*[, the HL equation (|l.lj) reads 

d t p = -b(t)d a p; 
P > 0; 

p(o,-) = r> ; 

{ D(p(t)) = . 
The above system reduces to 

p(t,a) = p ((T-x(t)) ; U7) 

D(p(t)) = 0, for all t in [0;t»] . 1 ' ' 
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The second equation in Q4.2|) is compatible with the first one as long as 

/ Po = , for all t in [0; i#] . 

J\a+x(t)\>l 

Therefore there exists a maximal time interval [0; T c ] on which the HL equation may reduce 
to a mere transport equation and this is for an intrinsic time T c (possibly infinite) defined by 



T c = M{t>0; p >0.} . (4.3) 

[ J\*+x(t)\>i J 

Note that T c is completely determined by the data po and b. If T c = +oo, the steady state 
pit, a) = po(cr — x{t)) i s a solution of the HL equation for all time. We shall now exhibit 
circumstances under which it is not the unique solution. For convenience, we restrict ourselves 
to the case when 6 = (we then have obviously T c = +oo). 

For po G L 1 (R) n L°°(R) such that p > 0, let us denote by F po the function from R+ to R + 
defined by F po (0) = D(po) and by 

Vx > 0, F po (x) = a / ( / p (a')ip/^-(a - a') da' ) da . 

J\a\>i \Jn J 



Proposition 4.1 Let po satisfy (|1.5[) and be such that -D(po) = 0, then 
i. If F Po satisfies 



J 

Jo 



+oo, (4.4) 



then p(t, a) = Po(a) is the unique solution to (|1.1|) in C 4 (L^) ; 

ii. Otherwise, (jl.lj) has an infinite number of solutions in C^(L^). The set of solutions to 
(|1.1|) is made of the steady state p(t,a) = poia) and of the functions (qt )t >o defined 
by 

Po(o-) ift<t 
q(t-t ,a) ift>t 



Qt (t,a) 



where q is the unique solution to (|1.1[) in C t °(L^) such that D(q) > on ]0, +oo[. 
Besides, 

p € — *q strongly in Lf loc (Ll). (4.5) 



e^0 



Lemma 4.1 Let p £ L X (1R) n L°°(R) such that 

Po>0, [ Po = l, D{p ) = 0. 

The function Fq is in C°([0, +oo[) D C°°(]0, +oo[), and is positive on ]0,+oo[. In addition, 
F po > on ]0,+oo[. 
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Proof of Lemma gUJ It is easy to check that F po £ C°([0, +oo[) n C°°(]0, +oo[), and that 
F Po > on ]0, +oo[. Since D(po) = 0, the function po is supported in [—1,1]. Thus, for any 

x > 



n / Po[<r) / n i- r 



-da da' 



l 

a ^ p (a' 



-l+cr' fi -cr 2 /4x /-+00 e -cr 2 /4x 



_ da + . _ _ 



da da' 



a 



+oo 



e"* 2 di + / e~* 2 di ) dcr' . 



\-<T< 



It follows that for any x > 0, 



1 



^ (x)=a ^y-/ o(fj/) v 4,3/2 



; . /l+cr' dW) 2 1-0-' d-^) s 



g 4se +■ 



4 x 3/2 



■e 4x 



da' > . 



(4.6) 







Lemma 4.2 Let 7 > and p G L^R) n L°°(R) such that 

Po>0, p = l, / \a\p < +00, £>(po) = 0. 

Jr Jn 

Let us consider the problem 

dtw = D(w(t))d^ a w — 7U> 
w(0,a) =p (a). 



(4.7) 



i. If F po satisfies \4-4\) then p(t, a) = Po(a) is the unique solution to 14.7]) in C^(L^); 

ii. Otherwise, A4-7\) has an infinite number of solutions in Cf(L^). The set of solutions to 
l4-7j) is made of the steady state w(t,a) = Po(a) and of the functions (vt ) to>0 defined 
by 

j v(t-t ,a) ift>t 
where v is the unique solution to (|4.7jl in C^(L^) such that D(v) > on ]0, +oo[. 



Corollary 4.1 The initial data po = 5X]-l,l[ fulfills the assumptions of the above lemma and 
io F dX (x) ^ +°°- Therefore there are infinitely many solutions to the equation (|1.4|) in the 



1 P0 

introduction. 
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Proof of Corollary I4.lt The only point to be checked is that F dx / X \ < +00. With the 



standard notation erfc(-z) = J z +0 ° e <2 dt, and by using q4.6|) and symmetry considerations, 
simple calculations yield 



Fp (x) 



1 

2a^x f 



TT 



erfc(d) da 







2a 



Since erfc(z) ~ ^ e z / z fox z going to +00, F po (x) ~ near and the integrability of 



l/F po on [0;1] follows. 



<> 



Proof of Lemma 14.21 

Let us consider a non-neg ative function D £ C°([0, +oo[, R + ). The unique solution in C^(L^) 



of the problem 



is given by 



d t w D = D{t)dl a w D - >yw D ; 
w D (0,a) =p (a), 



(4- 



w D (t,a) 



e 7 *po(^) 



7 * / PoOOv 



if t < t* 

(a - a') da' if t > t* 



(4.9) 



where t* = inf |i > 0, J D > Oj. Any solution to (|4.7|) thus satisfies 10 = w D ^ and 



therefore 



a 



ae 



H>1 

7* 



WD(w){t,o) da 



\a\>l \Jn 



-(a — a') da' ) da 
) 



y/2j* D(w(s)) da 

= e-^F po ^D(w(s))d s y 
It follows that the function D(w) is solution in C°([0, +00 [) to the nonlinear integral equation 

y(t)=e-^F po ^y(s)d s y (4.10) 

On the other hand, if D 6 C°([0, +oo[) is solution to (|4.10|) it is easy to check that the 
function wd defined by Q4.9JI is solution to (|4,8|) . 

If condition ()4.4|) is fulfilled, equation (|4.1(Jf) has a unique solution in C°([0, +00 [) (the constant 
function equal to zero) and the steady state w(t,-) = po thus is the unique solution to (|4.7|) 
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in C t °(L 2 ); otherwise, the set of solutions to (|4.1Uj) is made of the steady state w(t,-) = po 
and of the family (yt )t >o with 



Vto(t) = 

where the function z is defined on [0, +00 [ by 



if t < t 

z(t - t ) if t > t 



z{t) dx j if 7 >0; 

F(x) \ t , otherwise . 

Statement ii. is obtained by denoting by v the solution to (|4.8|) associated with the function 
z(t). 

Proof of Proposition I4.lt 

The solution p e to equation (|1,6|) satisfies the inequalities 

p~(t,a) <p e (t,(j) <pf(t,a) almost everywhere 
where p~ and pf are defined in C 4 °(L 2 ) by 

dtp; = (D(p e (t)) + e) dl a p~ e -p-; \ d t pt = (D( Pe (t)) + e) d? a p+ + ^-5 ; 

p 7 (0,a)= P0 (a); \ p 7 (0,a) = p (a). 

Therefore on the one hand 

D(p € (t)) > D(p-(t)) = e- l F po Qf (D(p 6 ) + e )) (4.11) 

and on the other hand 



< F„ Qr (%) + e)~) + i £ D(r)(») is. 



If (|4.4I) is not fulfilled, using (|4.11|) and the property that F po is strictly increasing on [0, +00 [, 
we obtain that 

D( Pe ) > Z(t) 

where z(t) is the function defined in the proof of Lemma 14.21 As for any < to < T, there 
exists 77 > such that z(t) > rj on [to,T] the same reasoning as in the non-degenerate case 
leads to the conclusion that (p e ) converges up to an extraction to p in T>'(]0, +oo[xE) and in 
L 2 ([t ,T],L 2 (R)) for any < t < T < +00, p being a solution to (jUJ) in C°(]0, +00 [, l? a ) 
such that D(p) > on ]0, +oo[. 
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5 Steady states 



Throughout this section the shear rate 6 is assumed to be a given constant and we are looking 
for solutions in L 1 (R) to the following system : 



D(v) 

-b d a p + D{p) d 2 aaP - XR\[-i,i] V + Sq((t) = on (0; T)xR 



P > 
D(p) = 



a 



p = 1 ; 
R 

a / da 

J\a\>l 



(5.1) 
(5.2) 

(5.3) 



Our main results are summarized in the following : 
Proposition 5.1 

i 7/6 = 0, any probability density which is compactly supported in [— 1;+1] is a solution 
to (|5.1[) which satisfies D(p) = 0. If a < |, i/iese are i/te on/y stationary solutions 
(and there are infinitely many), whereas when a > ^ £/iere exists a unique stationary 
solution corresponding to a positive value of D, which is explicitly given by (|5,4j) and 
(j5.fi)) below. This solution is even and with exponential decay at infinity. 

ii If b ^ 0, for any a > 0, there exists a unique stationary solution to 1)5.1)1 . and it 
corresponds to a positive value for D, which is implicitly given by ()5.7)) and ()5.8)) below. 
This solution has exponential decay at infinity. 

Remark 5.1 The statement in the above proposition is already pointed out by Hebraud and 
Lequeux J^. 



Proof of Proposition 15. H 
The case when 6 = 

We first observe that any non-negative function p which is normalized in L 1 (M) and with 
support in [—1; +1] is a solution to the system 1)5.1)1 since in that case all terms in equation 1)5.1)1 
cancel. We now examine the issue of existence of solutions of 1)5.1)) such that D(p) > 0. For 
simplicity we denote D = D(p). For given constant D > 0, it is very easy to calculate 
explicitly the solutions of 1)5.1)1 on each of the three regions a < —1, a G [— 1; +1] and a > 1. 
Using compatibility conditions on R and the fact that p has to be in L 1 (R) one obtains : 



if a < -1 , 
if - 1 < o- < , 
if < a < 1 , 
if 1< a . 



p(a) 





f<r)/-/D 


2a 




1 l v 


^D + l 


2a a + 


2a 




1 i_ v 


/ Z) + 1 




2a 








2a 





(5.4) 
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The compatibility condition D = D(p) happens to be then automatically satisfied and the 
normalization constraint j-^p = 1 imposes that D solves 

D + V15 = a-^. (5.5) 

Since D > 0, we immediately reach a contradiction when a < g, whereas when a > \ 
equation Q5.5|) admits a unique positive solution; namely 



D 



1 



+ 



Ma - 1 



(5.6) 



2 2 
The case when 6^0 

First of all, we observe that if D = every term in equation 1)5.1)1 but vanish. Thus 

p has to be a non-zero constant which is in contradiction with p £ L X (1R). So necessarily 
D > 0. For given positive constant D, we then solve 1)5.1)1 as above and obtain 

,(3+a 



Pip) 



a\ e' 

Jl fj D 
a 2 eo + a 2 - — 



«2 



ai e 



+ a 2 



if a < -1 , 
if - 1 < a < , 

if < a < 1 , 
if 1 < a , 



(5.7) 



with 



8 ± = ± - 

H 2D ' 



1 b 2 + AD 



D 2 



' D 



and 



ai 



CL2 



a 



e b/2D _ (3 - e ~b/2D^ ) 

D(3 + e b / 2D 



ab (p+e b / 2D - p-e~ b / 2D ) ' 

It is tedious but easy to check that this function always fulfills the self-consistency condition 
D = D{p) and that the normalization condition J^p = 1 reads 

D {l + p + ) + (f3- -l)e - b l D 
T B+ - 8- er b / D 



+ D 



a 



/3+ _ (3 - e - b/ D ' -■ ( 5 - 8 ) 

For any b > (the negative values of b are dealt with by replacing a by — a), the left-hand 
side of (|5.8[) is a continuous function which goes to +oo when D goes to infinity and goes to 
zero when D goes to 0. This already ensures the existence of at least one steady state for any 
a > 0. Moreover, setting z = (for example) we may rewrite the left-hand side of 1)5. 8|) as 

1 + ^coth(z/26) + ^(z 2 + 4z)V2" 



/(*) 



62 2&2 

z z 



z+(z 2 + 4z)V2 coth(z/26) 



Next we check that the function / is monotone decreasing (thus, the left-hand side of 1)5.8)1 
is increasing with respect to D), whence the uniqueness result. (} 
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